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A potential flux of high-energy neutrinos from the annihilation of dark matter particles trapped 
within the Sun has been exploited to place indirect limits on particle dark matter. In most models, 
the dark matter interacts weakly, but the possibility of a dark matter particle with a large cross 
section for elastic scattering on other dark matter particles has been proposed in several contexts. 
I study the consequences of such dark matter self-interactions for the high-energy neutrino flux 
from annihilation within the Sun. The self-interaction among dark matter particles may allow dark 
matter in the halo to be captured within the Sun by scattering off of dark matter particles that 
have already been captured within the Sun. This effect is not negligible in acceptable and accessible 
regions of parameter space. Enhancements in the predicted high-energy neutrino flux from the Sun 
of tens to hundreds of percent can be realized in broad regions of parameter space. Enhancements 
as large as factors of several hundred may be realized in extreme regions of the viable parameter 
space. Large enhancements require the dark matter annihilation cross section to be relatively small, 
((tad) < 10 _27 cm 3 s _1 . This phenomenology is interesting. First, self-capture is negligible for the 
Earth, so dark matter self-interactions break the correspondence between the solar and terrestrial 
neutrino signals. Likewise, the correspondence between indirect and direct detection limits on 
scattering cross sections on nuclei is broken by the self-interaction. These broken correspondences 
may evince strong dark matter self-interactions. In some cases, self-capture can lead to observable 
indirect signals in regions of parameter space where limits from direct detection experiments would 
indicate that no such signal should be observable. 

PACS numbers: 95.35.+d,95.30.Cq,95.55.Vj,98.35.Gi,98.80.Cq 



I. INTRODUCTION 

A great deal of observational evidence indicates that 
a form of non-relativistic, non-baryonic matter consti- 
tutes the vast majority of mass in the Universe. The 
unknown nature of the dark matter that binds galaxies 
and drives cosmic structure formation remains an impor- 
tant problem in cosmology and particle physics. Among 
dark matter candidates, weakly-interacting massive par- 
ticles (WIMPs), including the lightest superpartners in 
supersymmetric theories, have received the most atten- 
tion (for a review, see Ref. In this paper, I study 
a potential enhancement in high-energy neutrino fluxes 
from dark matter annihilations within the Sun in models 
where a WIMP-like dark matter particle exhibits rela- 
tively strong interactions with itself. 

Indirect, astrophysical probes of dark matter are an 
important element of any comprehensive program to 
identify the dark matter unambiguously. One indirect 
probe of WIMP dark matter is a potentially-detectable 
flux of high-energy muon neutrinos arising from the anni- 
hilation of dark matter particles captured within the Sun 
@, 1 1 IE S, 0, S, i M El, El. A similar signal from 
within the Earth may also be exploited in this regard 
[E IE EI, and though the terrestrial signal is typically 
smaller than the solar signal, it is a valuable cross-check 
[Io| . In fact, these signals have already been brought to 
bear to limit dark matter elastic scattering cross sections 
with nucleons at interesting levels [3, ELEE EE Ell • 

This basic scenario is simple. As the Sun moves 
through the halo of WIMPs, some of the WIMPs scatter 
elastically off of nuclei within the Sun. Many WIMP- 



nucleus interactions result in WIMPs moving at speeds 
lower than the local escape speed relative to the Sun. 
These particles are captured and for a large region of rel- 
evant parameter space they come to thermal equilibrium 
in the interior of the Sun. Eventually, the build-up of 
WIMPs within the Sun is limited by the annihilation of 
these WIMPs producing neutrinos that can escape from 
the Sun. Annihilation products other than neutrinos in- 
teract within the Sun and are not observable at Earth. 

Dark matter particles that interact weakly with stan- 
dard model particles, but exhibit comparably rather 
strong interactions among themselves have now been pro- 
posed in several different contexts [lj5|, HE HE HE El, HE 
M, HE HE pi HE HE HE HI. Some bounds on dark 
matter self interactions exist [33t HE HE HE H3] and ob- 
servational tests of various scenarios have been proposed 
[H, HE HE EL 52 , but a wide range of parameter space 
remains and will remain viable for M x ~ a few x 10 2 GeV 
dark matter particles with large self-interaction cross sec- 
tions, er xx ~ 10~ 24 cm 2 . 

Large cross sections for dark matter particles to scatter 
elastically off of each other open a new possibility for 
capture within the Sun. In addition to nuclei, previously- 
captured dark matter particles or dark matter particles 
otherwise sequestered within the solar interior may serve 
as additional targets for the capture of halo dark matter 
particles. I refer to this as dark matter "self-capture" and 
it is this possibility that I consider in detail in this paper. 
Previous studies have considered distinct modifications 
to high-energy neutrino fluxes from the Earth and Sun 
due to inelastic scattering of dark matter with against 
nuclei [H,[y|. 
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I begin in § [II] with a brief sketch of the standard 
scenario for indirect detection of dark matter via high- 
energy neutrinos from the solar interior. In § Mil I use 
a simple, order-of-magnitude estimate to show that dark 
matter self-capture within the Sun may not be negligi- 
ble for a range of viable and interesting models. On the 
other hand, self-capture within the Earth is always negli- 
gible for models that have not yet been excluded by other 
means. 

In § IIV1 I give the results of more detailed calculations 
of the importance of dark matter particle self-capture 
within the Sun. These results demonstrate that mod- 
est enhancements of tens to a hundred percent relative 
to models in which self-capture is negligible are possi- 
ble over reasonably broad ranges of interesting parame- 
ter space. Significantly larger flux enhancements of up 
to factors of hundreds are possible in extreme regions of 
the dark matter parameter space. Throughout, I remain 
relatively agnostic about the nature of the dark matter 
and present results as a function of the most directly 
relevant model parameters, dark matter particle mass 
M x , dark matter-proton scattering cross section cr p , dark 
matter self-interaction cross section cr xx , and thermally- 
averaged dark matter annihilation cross section multi- 
plied by relative velocity (<jav). However, I do use the 
findings of detailed explorations of the parameters avail- 
able to neutralino dark matter in supersymmetric scenar- 
ios as g uidance for interesting values of these parameters 

I summarize my results and conclusions in § [V] In 
particular, I emphasize that flux enhancements from self- 
capture scenarios may be important for two reasons. 
First, the solar flux may be significantly altered by self- 
capture, while the terrestrial flux cannot be. Therefore, 
the ratio of the solar to terrestrial high-energy neutrino 
fluxes from dark matter may be markedly different from 
the standard predictions and this may signify new dark 
matter interactions. Likewise, the similar correspon- 
dence between direct search signals and solar high-energy 
neutrino fluxes can be broken. In some cases, models 
that may otherwise be ruled out by direct dark matter 
searches may produce observable neutrino signals due to 
the self-capture enhancement. I include in an Appendix 
the details of the capture rate calculations that I perform, 
including an example of the capture rates that I compute 
in the standard scenario of spin-independent capture off 
of nuclei. This discussion follows the derivations given 
by Gould [HI El- 



II. HIGH-ENERGY NEUTRINOS FROM THE 
SUN AND DARK MATTER SELF-CAPTURE 

In the most well-studied scenarios, captured dark mat- 
ter particles typically thermalize in the solar interior 
on a timescale less than the age of the Sun (r Q w 
5 x 10 9 yr) as well as the other timescales in the problem 
@, i, H i, i, 0, H El- In this case, the time evolution 



of the number of dark matter particles in the Sun N x , 
follows 



dN x 
dt 



C c + C S N X - C a N> 



(1) 



The coefficients are the rate of capture of dark matter 
particles by scattering off of nuclei within the Sun C c , 
twice the rate of annihilation per pair of dark matter 
particles within the Sun C a {twice because each annihi- 
lation eliminates two particles) , and the rate of capture 
of dark matter particles by scattering off of other dark 
matter particles that have already been captured within 
the Sun C s . The C S N K term on the right-hand side of 
Eq. (JT|) is the new term that I study in this paper. Cap- 
ture rates were first computed by Press and Spergel .9] 
and this calculation was revised, corrected, and greatly 
expanded upon in an impressive series of papers by Gould 
[HI EH El • In principle, evaporation of captured parti- 
cles from the Sun can also occur, but this is unimportant 
for masses larger than a few GeV 0, EH ■ I discuss the 
specific rates at greater length below and in the appendix. 
For the time being, let us focus attention on Eq. |T|). 

In the standard treatment, self-capture of dark matter 
particles is ignored (C s = 0). The solution of Eq. ([1]) for 
N x = at t = is then 



tanh(v/C c C7 a i). 



(2) 



There is a timescale for equilibration between dark 
matter annihilation and dark matter capture, r oq = 
l/y/C c C a . For many models of interest, r cq << r© and 
the solution approaches a steady state, 



(3) 



In this circumstance, the rate of annihilation of captured 
dark matter particles within the Sun is 




(4) 



The factor of 1/2 in Eq. (H)) arises because there are 
Nl e j2 distinct pairs of particles, while this factor is not 
present in Eq. {T]) because each annihilation eliminates 
two dark matter particles from the Sun. The annihila- 
tion rate T a is independent of the annihilation rate coef- 
ficient C a . Consequently, the observable flux at Earth is 
independent of the dark matter mutual annihilation cross 
section, provided the cross section is large enough that 
the equilibrium solution [Eq. ||3J)] obtains. Dark mat- 
ter particles annihilate upon capture and the flux at the 
Earth is modulated only by the capture rate C c . 

In the first papers to study high-energy neutrinos from 
dark matter annihilation within the Sun, it was noted 
that this phenomenology is interesting and useful. First, 
the flux from annihilations at rate Eq. ^ is independent 
of annihilation cross section, so this is an indirect search 
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FIG. 1: The evolution of the annihilation rate of dark mat- 
ter particles captured within the sun as a function of time. 
The solid lines shows evolution to the steady-state solution of 
Eq. (|3} and Eq. in the absence of significant dark matter 
self interactions. The dashed lines show evolution to the new 
equilibrium in a model in which self-interactions are impor- 
Time is shown in units of the stan- 
and the squared 



tant and C s = 

dard equilibration timescale r cq = 
number of captured dark matter particles is shown relative to 
the standard equilibrium number iV x , eq = yCc/CsL- 



method that does not rely on models with relatively large 
annihilation cross sections. Second, the high-energy neu- 
trino signal from annihilation in the solar interior can 
be related to other potentially-observable signals. I dis- 
cuss the rate of capture in the appendix, but it suffices 
to note that the capture rate of dark matter particles in 
the Sun should be proportional to the cross section for 
a dark matter particle to scatter off of a nucleus (on) 
within the Sun and the local density of dark matter (p x ), 
C c oc ctnPx- A similar process may operate within the 
Earth, whereby captured dark matter particles give rise 
to high-energy neutrinos from the Earth's center and this 
signal should also grow in proportion to the product ctn/?x 
@H, Ell- Moreover, signals in direct dark matter search 
experiments are proportional to this same product as well 
Q. As a consequence, indirect detection of dark matter 
through high-energy neutrinos from the Sun and Earth as 
well as direct dark matter search experiments may serve 
to check each other and corroborate any detections. 

If self-capture is not ignored this picture may be al- 
tered. In particular, the general solution to Eq. (JTJ) with 
C s > is 



C c tanh(t/C) 



where 



C- 1 - C s tanh(i/C)/2 



(5) 



(6) 



The steady-state solution at t ^> C is 




(7) 



Clearly, Eq. Q reduces to Eq. in the case that C s = 0. 
Whether or not self-capture of dark matter particles is 
important can be discerned by comparing the timescales 
C~ l and 1/VCcCa, and I will pursue this comparison 
shortly. Consider the case where Cg 3> C c C a . In this cir- 
cumstance, iV^ e eq ~ C s /C a . The annihilation rate within 
the sun is then 



1 



C 2 JC & . 



(8) 



The annihilation rate grows in inverse proportion to the 
rate coefficient C a when self-capture is possible. As I 
discuss below, the annihilation rate coefficient is pro- 
portional to the annihilation cross section, so models 
with relatively low annihilation cross sections are favored 
for an indirect neutrino signal from the Sun when self- 
capture is not negligible. 

Figure [1] shows two examples of the evolution of the 
rate of annihilation of dark matter particles captured 
within the Sun as a function of time. If self-capture is 
not negligible, the evolution of the total number of cap- 
tured dark matter particles may come to be dominated 
by the second term on the right-hand side of Eq. |1} and 
N x may grow exponentially for some time. This expo- 
nential growth is eventually terminated by annihilations, 
so the lower the annihilation cross section, the more ex- 
ponential growth of N x is important and the greater the 
relative flux of neutrinos from the Sun may be (provided 
the Sun is old enough that the equilibrium solution has 
been achieved). 



III. THE IMPORTANCE OF DARK MATTER 
SELF-CAPTURE 

I now present a simple estimate of the parameter val- 
ues for which dark matter particle self-capture may be a 
non-negligible effect. I define the ratio i? s = C^/C c C a . 
Comparing to Eq. ((SJ) and Eq. (J7J), self-capture will be 
negligible when R a <C 1 and dominant when R s 3> 1. 
The task now is to evaluate each of the rate coefficients 
for a set of parameters describing the interactions of the 
dark matter particle. 

I show in the Appendix that a simple approximation 
for the self-capture rate is 



C S N X 



2 ^x^xx^csc 



(R©) 



-iVx(^x) 



rf (77) 
V 



(9) 

where n x is the local number density of dark matter par- 
ticles in the halo, <7 XX is the dark matter elastic scat- 
tering cross section, v CS c(R-o) is the escape speed from 
the surface of the Sun, v is the local three-dimensional 
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velocity dispersion of dark matter particles in the halo, 
(</> x ) ~ 5.1 is a dimensionless average solar potential ex- 
perienced by captured dark matter particles within the 
Sun, and rj 2 — 3(vq/v) 2 /2 is the square of a dimension- 
less speed of the Sun through the Galactic halo. Relation 
Eq. (J9]) neglects the recoils of the target dark matter par- 
ticles; however, I demonstrate in the Appendix that this 
is a reasonable approximation when escape speeds are 
significant compared to vq and v. This condition holds 
for capture within the Sun. 

The full expressions for capture off of nuclei are un- 
wieldy, so I make an effort here to evaluate R s ap- 
proximately. To simplify the evaluation of R s , I tem- 
porarily assume that capture of dark matter by scat- 
tering off of nuclei is never significantly restricted by 
the kinematics of the scattering. This is not gener- 
ally true. In any individual scattering event dark mat- 
ter particles can lose a fraction of their kinetic energy 
AE/E < 4M X M N /(M X + M N ) 2 , where M N is the mass 
of the target nucleus. Capture within the Sun typically 
requires AE/E > 1/5. Therefore, capture can be effi- 
cient for dark matter and nucleus masses that differ by a 
factor of as much as ~ 20. If the mass of the dark matter 
particle is sufficiently different from the mass of the nu- 
cleus on which it scatters, capture may be kinematically 
unfavorable. 

Neglecting kinematic limitations to capture on nuclei, 
the capture rate off of a particular nuclear species N is 
given by a relation analogous to Eq. |J9]) with C s iV x — * C c , 
ct xx -> cr N , N * -* /nM /M n , and (<j> x ) -> (0 N ). The 
quantity ctn is the scattering cross section of the dark 
matter particle off of the nucleus N, /n is the fraction 
of the solar mass in nucleus N, and (0n) is the average 
dimensionless potential experienced by these nuclei. For 
most nuclei within the Sun, (4>n) — 3.2 [10j]. Examining 
Eq. ([5]) , the capture rate off of nuclei in this limit scales 
with dark matter mass as C c oc M~ l . When M x ^> 
Mn, the kinematic limitation to the capture rate is not 
negligible, the maximum fractional kinetic energy lost 
per collision is ~ 4Mn/M x , and the capture rate scales 
as C c oc M~ 2 . 

Assuming that the dark matter particles equilibrate 
with the solar interior rapidly upon capture, the coeffi- 
cient C a is likewise simple to estimate. Let e x (r) be the 
number density of dark matter particles as a function of 
position in the Sun. The annihilation rate coefficient is 
then 

a = ^f% x M^, (io) 

and the naive expectation for a non-relativistic thermal 



relic dark matter particle is that (<tav) ~ 10~ 26 cm 3 s _1 . 
Under the assumption of a thermal distribution at an 
effective solar core temperature T . c = 1.57 x 10 7 K, 
the distribution e x (r) oc exjp[— M x (j)(r) /Tq_ c ], where </>(r) 
is the gravitational potential as a function of position 
within the Sun. Making a further assumption of a con- 
stant solar density of p© iC = 150 g/cm 3 the integral 
is straightforward to evaluate. Conventionally, this has 
been represented in terms of effective volumes 

C & = {a A v)^ (11) 

where 

^^TSxlO-^^y^cm 3 . (12) 

The effective volumes represent ~ 10~ 6 of the total solar 
volume, so captured particles extend over only ~ 1% of 
the solar radius, justifying the constant density approxi- 
mation [jj. 

The question arises whether the strong self-interactions 
among the dark matter particles should alter the assump- 
tion of a thermal distribution. A definitive answer to this 
question requires solving the Boltzmann equation. Such 
a calculation is extensive and beyond the scope of this 
paper. It is reasonable to suspect that any modifications 
will be minor for parameters of interest. In the relevant 
parameter regime (which will be more clearly delineated 
below) the increase in the number of captured dark mat- 
ter particles relative to the standard scenario is modest 
(less than a factor of ~ 30), yet collisions among dark 
matter particles may happen at a rate that is compara- 
ble to the rate of collisions between dark matter particles 
and nuclei. It is a relatively s imp le matter to use the ap- 
proximate methods of Ref. [12j to show that particles 
with M x > 10 GeV remain localized well within the solar 
interior and are not altered by the temperature gradient 
within the Sun, and that the rate of energy inflow due 
to capture is slower than the rate of thermalization with 
the solar interior. Finally, it is also straightforward to 
show that for masses < 10 12 GeV (which far exceeds the 
unitarity bound for a thermal relic, e.g. [48]), the dark 
matter density is never sufficiently large for the dark mat- 
ter to be self-gravitating, so the self-interaction does not 
lead to rapid collapse via the gravothermal catastrophe 
(the case of M x > 10 12 GeV is treated in Ref. [H|). 
Given these considerations, a standard thermal profile 
seems a reasonable approximation and I proceed under 
this assumption. 



With expressions for C c , C s , [Eq. ©] and C a [Eq. (fTTj) ], all of the pieces are now in place to approximate the ratio 
R s = C 2 /C c C a . Making the simplifying assumption that capture occurs primarily off of a single type of nucleus, this 
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gives 

R _ /3 cr 2 x ?; csc (R ) i> csc (Rq) (0 X ) 2 Vi_ n x M^Vi erf(?y) 
s ~ V 2 a N {a A v} v (0 N ) y 2 / N M n 

Evaluating this for the particular case of capture off of Oxygen (which is the most important individual element for 
dark matter capture in the Sun, see Figure [5] in the Appendix and Ref. [13]), taking /n = 10 -2 (almost twice the 
solar Oxygen abundance to account for the simplicity of the current estimate) , and keeping relevant aspects of the 
dark matter particle model explicit yields, 



The result in Eq. 1|14|) is somewhat startling. Cur- 
rent best bounds on the elastic scattering cross section 
of dark matter particles off of each other are assumption- 
dependent and approximate, but they indicate that 
o-xx < 10" 23 (Afx/100 GeV) cm 2 [H l| H H, H| . A 
slightly less restrictive bound from an analysis of the bul- 
let cluster is probably the least dependent upon particu- 
lar assumptions [33| . Eq. (TTJ} indicates that dark mat- 
ter particle self-capture within the Sun is not necessarily 
a negligible effect in acceptable regions of dark matter 
particle parameter space. Moreover, Eq. (fT4| neglects 
the fact that dark matter capture off of nuclei may be 
kinematically unfavorable, while self-capture can never 
be kinematically unfavorable because the particles will 
always have equal mass. 

I note in passing that self-capture by the Earth can 
never be important. As 1 discuss in the Appendix, col- 
lisions between halo dark matter particles and particles 
already captured within a body may result in the target 
dark matter particles being ejected from the body upon 
recoil. The net result in this case is no gain in the num- 
ber of captured dark matter particles. Ejection by recoil 
depends upon the ratio of the speed of the particles at 
infinity to the escape speed from the body. Within the 
Sun, escape speeds are always significantly higher than 
the typical relative speeds of dark matter particles at in- 
finity and ejection is only a small correction to the simple 
solar capture estimate. In the case of the Earth, escape 
speeds are more than an order of magnitude lower than 
the typical relative speeds of dark matter particles at in- 
ifinity so almost all collisions within the Earth result in 
ejection of the target dark matter particle. The ejection 
of the targets from the Earth introduces the possibility 
that the halo dark matter particles may scour the Earth 
of any particles captured through interactions with nu- 
clei. This can be computed in a manner analogous to 
self-capture, though the sign of the term linear in N x in 
Eq. (fTJ) would be negative. In the case of the Earth, the 
ejection rate is small and R s < 10 -3 for all parameters of 
interest. Modifications to the Earth signal are negligible. 

In the following section, I show results from a more de- 
tailed calculation of the importance of dark matter parti- 



cle self-capture for high-energy neutrino fluxes observed 
at the Earth. I use the formulae from Ref. [10| to compute 
the capture rate of dark matter particles from nuclei as 
described in the Appendix. These formulae are lengthy 
and I do not reproduce them in full here, though I give 
an example of the capture rates that I use in Fig. [5] I use 
Eq. (|A.20|) to compute the rate coefficient for self-capture 
of dark matter particles. This relation is derived in the 
Appendix and includes the reduction in the capture rate 
due to the potential ejection of target dark matter parti- 
cles. 

The most relevant quantity to compute is the enhance- 
ment in the neutrino signal due to self-capture of dark 
matter relative to the neutrino flux expected in the stan- 
dard calculation. I define the quantity 



N 2 



iV 2 



(15) 



which is the ratio of the high-energy neutrino flux when 
self-capture is possible to the high-energy neutrino flux 
without the possibility of self-capture. My primary re- 
sults are illustrations of the dependence of (3 on the 
parameters er X x, (o~av), the spin-independent dark mat- 
ter particle-proton elastic scatting cross section cp 51 , and 
M x . To evaluate f3, I do not assume that the equilib- 
rium solutions of Eq. (fTJ) are attained. Rather, I evaluate 
N x (t — tq) and N KjS (t = Tq) using the general solution 
ofEq. ©. 

I assume that the cross section for dark matter scat- 
tering off of nuclei other than Hydrogen is given by [l[ 



M 2 M 2 



(M x + mproton) 2 



(M x + M N ) 2 M 2 m 



(16) 



proton 



where A is the atomic mass number and m pro ton is the 
proton mass. Loss of coherence is accounted for in the full 
formulae through suppression by an exponential form fac- 
tor [ToL fill] . In the following examples, I focus primarily 
on spin-independent interactions. Spin-dependent cap- 
ture of dark matter off of nuclei occurs only for Hydrogen 
within the Sun and is typically down by roughly two or 
more orders of magnitude at fixed cross section for high- 
mass (M x > 100 GeV) dark matter particles. This may 
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be mitigated by the fact that the spin-dependent cross 
section for scattering off of protons is typically ~ 1 — 2 
orders of magnitude larger than the spin-independent 
cross section in viable regions of the constrained min- 
imal supersymmetric standard model parameter space 
[ill ]. Moreover, direct search experiments typically use 
large nuclei with no net spin and exploit the scaling of 
Eq. (flB)) , so neutrino telescopes [3, [H, UH a very com- 
petitive with direct search bounds on a spin-dependent 
interaction [El, [EH, [El, [EH and should remain so [EH- 
In § IIV1 I show estimates for the spin-independent case 
as it is more general, including capture off of all nuclei 
within the Sun, and more interesting for present pur- 
poses in the sense that complementary constraints from 
direct search experiments are more competitive with in- 
direct methods for spin-independent capture. Including 
spin-dependent capture would typically add a term that 
is at most comparable to C c (though this is a model- 
dependent statement) and I find comparable values of 
f3 for spin-dependent capture with spin-dependent cross 
sections a^ D ~ lOVp 31 . 

It is important to set the scale of the signal relative 
to current and future observations, so I present esti- 
mates of absolute fluxes in § IIVI as well. Annihilations 
in the Sun, lead to a flux of high-energy neutrinos at the 
Earth. The observable signal at a detector such as Ice- 
Cube [HI, [El, [El| is a flux of upward-directed high-energy 
muons induced by scattering of muon neutrinos near the 
detector. The muon flux at the detector is therefore a 
relatively complicated product and may be written as 

* - sS^ / dE » f iE » 

x f E " d E,P(E^E»,\)^jM 

Je„ OMjfx 

x £ P osc (i^ M ) £ B[ ^l. (17) 

i £ 

Individually, the factors in Eq. (|17p are relatively simple 
to interpret. T a is the annihilation rate in the solar inte- 
rior, A§ is the semi-major axis of the Earth's orbit about 
the Sun, and «t is the number density of target nuclei 
near the detector. The quantity P(E^ — ► E^,X) is the 
probability for a muon of initial energy to have final 
energy E^ after traversing a path of length A in the detec- 
tor material, the differential cross section da u ^(E u )/dE^ 
describes the production of a muon of initial energy E^ 
from an incident neutrino of energy E y , P sc(i - * A*) is 
the probability that a neutrino produced as flavor i is 
a muon neutrino near the detector, Bf is the branching 
ratio to annihilation channel /, and dN^f/dE^ is the dif- 
ferential spectrum of neutrinos of flavor i per unit energy 
dE v produced per f-channel annihilation. 

I have evaluated Eq. (jTTJ) for an experiment such as Ice- 
Cube [H, [El, [El| using the results of the WimpSim Monte 
Carlo simulations [E3| as available through the DarkSusy 
package [El|. I choose -Eth = 1 GeV in accord with 



the common convention for reporting results from neu- 
trino telescopes. An instrument like IceCube observes 
events above tens of GeV, so sensitivities quoted rela- 
tive to -Eth = 1 GeV depend upon an assumed spec- 
trum. I show flux normalizations for two simple choices 
of branching fraction. I show results for annihilation to 
W + W~ gauge bosons only (B w +w- = 1) as a simple 
approximation of fluxes that may be produced from an- 
nihilation of a typical neutralino and a slightly more op- 
timistic case of annihilation to t + t~ only (B T + T - = 1). 
Annihilation to t + t~ yields about three times higher flux 
than annihilation to gauge bosons throu gh m ost of the 
relevant dark matter particle mass range [54 , [E3 • 

IV. RESULTS FOR HIGH-ENERGY NEUTRINO 
FLUX ENHANCEMENTS 

I summarize results on the relative importance of a 
contribution from dark matter particle self-capture to 
high-energy neutrino fluxes from the Sun in the contour 
plots of Figure [2] and Figure [31 There are four param- 
eters of most immediate interest, namely cr xx , M x , cr^ 1 , 
and (ctaw). Fig. [2] displays contours of constant j3 in the 
A/ x -cr xx plane for several fixed values of tJp 1 and (<tav), 
while Fig. [3] shows contours of f3 in the M x -<7p plane for 
specific choices of <r xx and (<tav). 

Consider Fig. [21 which shows an interesting, gen- 
eral set of results. The shaded regions at the upper 
left represent values of ct xx that are already ruled out 
[l9l l33l . [H], [H, [11|; however, computing these bounds 
is complex and the position of the boundary in each 
case remains somewhat controversial [33l . |34| . |35| so this 
boundary should be regarded as approximate. Notice 
that boosts in neutrino fluxes of several tens up to 100% 
can be achieved for quite reasonable parameter values. 
Significantly larger boosts of up to (3 ~ 10 3 can be real- 
ized in extreme regions of the viable parameter space. 

Beyond this, several additional features of Fig. [2] are 
worthy of explicit note. Eq. (fl4)) indicates that lines of 
constant [3 on the M x -cr xx plane should run as er xx cx 

M x ■ In practice, the lines of constant (3 are somewhat 
more shallow than this because the approximate form of 
C c [Eq. used in the simple estimate of Eq. (TH| as- 
sumed favorable kinematics for scattering off of nuclei at 
all dark matter particle masses. In the case of favorable 
kinematics, C c cx M~ . However, as the dark matter 
particle and nucleus masses become less comparable, the 

capture rate tends to C c cx M~ 2 and lines of constant (3 

3/4 

become shallower, approaching cr xx cx M x . In addition, 
comparing the pair of panels (b) and (c) or (d) and (e) 
in Fig. [2] it is clear that the scaling cr 2 x cx <7n(<ta.u) from 
Eq. (TTJ} for fixed (3 at a particular M x is valid. This is 
sensible, because these cross sections serve only to scale 
the rates C s , C c , and C a (so long as the equilibrium so- 
lution is achieved). 

Also evident in Fig. [5] is that large enhancements may 
only be achieved when annihilation cross sections are 
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FIG. 2: Factors of flux enhancement in the M x -a xx plane. Each panel shows contours of constant relative flux enhancement 
in models of self-interacting dark matter. The panels (a)-(f) are labeled with assumed values of ap 51 and (<tav). The shaded 
regions at the upper left correspond to parameter values that are disfavored by analysis of either the Bullet Cluster [33|] or 
galaxy cluster shapes [3, l34l 13a . |36| ] . 



relatively low {crpjv) < 1CP 26 cm 3 s _1 , where the nu- 
merical value is the canonical value for a thermal relic 
dark matter particle. This can be seen most dramati- 
cally in panel (f) where I have taken er^ 1 = 10~ 45 cm 2 



and (<7av) = 1CP 26 cnr^s^ 1 and the enhancements are at 
most a few percent over viable parameter ranges. Even 
discounting experimental limitations, it is thought that 
the intrinsic errors in computing neutrino fluxes from 
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dark matter capture within the Sun should be a few 
tens of percent 0, HO, HI H M, [12, so this indi- 
cates that such an effect can only be interesting when 
(<tau) ^ a few x 10~ 27 cm 3 s _1 . Within the context of 
scans of restricted regions of supersymmetric parameter 
space, annihilation cross sections well below this value 
are achievable [45], H^, [47], [(53| , so significantly lower cross 
sections are attainable in theories with complex parti- 
cle spectra. Even in exceedingly simple proposals there 
exists sufficient freedom to set cr xx and (a^v} apart sig- 
nificantly [12, Consider interaction via exchange 
of a boson of mass my. The perturbative annihila- 
tion and scattering cross sections should be related as 
ca/cxx ~ (my/Mx) 4 and both my and the coupling 
strength remain to be fixed. However, these results do 
indicate that models of dark matter self-interaction that 
lead to very large annihilation cross sections (such as the 
Sommerfeld-enhanced scenarios of significant recent in- 
terest, see Refs. [H HH, [H, HI) will induce little addi- 
tional neutrino flux due to self-capture. 

Figure [3J shows contours of (3 in the M^-a^ plane for 
four choices of <7 XX and (ctav) and complements the re- 
sults in Fig. O First, the results of current direct dark 
matter searches can be compared in this plane. Di- 
rect search experiments constrain cp 51 directly. Current 
bounds place af > 10~ 43 cm 2 at M x - 10 2 GeV 
[Hot l5ll . l52l [53 | . This bound becomes slightly better with 
decreasing M x until M x ~ 50 GeV and at higher masses 
this bound grows ex M x . The competition is much less 
severe from spin-dependent searches [HI, [EH, [53| , as the 
indirect limit from the Sun already exceeds the direct 
search limit by more than an order of magnitude over 
a wide range of masses [l6|. Second, Fig. [3] shows con- 
tours of constant absolute muon flux above a threshold of 
^th = 1 GeV in models where there is no self-interaction 
in order to set the absolute scale. Two annihilation chan- 
nels are shown, annihilation into W + W~ (solid gray) and 
t + t~ (dashed gray). Current experiments are limited 
to muon fluxes above several hundred per km 2 per year 
[iH , 1551 [56j . Assuming the relatively hard spectra from 
annihilation to gauge bosons, IceCube with DeepCore 
extension should optimistically be capable of detecting 
fluxes down to ~ 60 km yr for particle masses above 
~ 200 GeV, with relatively lower sensitivity below this 
mass [13, [H, [H, [iH [65|. A deep-sea neutrino facility 
such as the KM3NeT effor t 161 . building on the ANT ARES 
[13, [Si, IH, [ZO, [ill , NEMO [Zalzllzl, and NESTOR [zllz^ 
work, may achieve comparable or better sensitivities. 

Figure [3J also shows regions where the equilibrium so- 
lution of Eq. ([7]) has not yet been attained for a solar 
age of tq — 5 x 10 9 Gyr. I approximate the equilib- 
rium boundary as the contour where the predicted flux 
is 58% of the value it would be at equilibrium, because 
tanh 2 (l) ~ 0.58 [see Eq. ([3]). Equilibrium is achieved in 
the majority of the paramete r sp ace corresponding to a 
potentially-detectable signal (54j. However, notice that 
the contours of constant flux at Earth in scenarios with 
no self-capture are not the same in each panel [partic- 



ularly so in panel (d)] because fluxes are no longer de- 
termined solely by up 51 and M x for models that are not 
equilibrated. 

Consider panels, (b)-(d) of Fig. [3J In these panels, the 
equilibrium boundary exhibits a very shallow minimum 
in <7p near M x of a few hundred GeV. In these panels, 
the equilibrium boundaries are essentially the same as 
they would be in the absence of any dark matter self 
interaction. The minimum occurs due to the compe- 
tition between capture and annihilation in the relevant 
timescale, r cq = l/\/C c C a . The annihilation rate scales 

3 /2 

with dark matter mass as C a cx M x [see Eq. (fTTT) and 
Eq. I|12p]. while at relatively low masses C c tx M~ x . For 
dark matter particle masses greater than several hundred 
GeV, the capture rate transitions to the regime where 
kinematic suppression of capture becomes important and 
C c cx M~ 2 (I have neglected the orbital effects that also 
tend to slow thermalization within the Sun for high-mass 
dark matter candidates [60|). 

On the other hand, panel (a) of Fig. [3J shows a distinct 
feature in the equilibration boundary at M x ~ 300 GeV. 
The feature is caused by self-capture. At high <r xx and 
low Mx, self-capture of dark matter is important and can 
drive rapid equilibration even for very low values of C c . 
In the absence of self-capture, the equilibration boundary 
in panel (a) of Fig. [3J would be relatively flat as a function 
of M x as in panels (b)-(d). 

The contours of constant (3 in Fig. [3J show that inter- 
esting regions of parameter space can lead to detectable 
boosts in muon fluxes at Earth of tens of percent to 100%. 
Somewhat more extreme choices of parameters can lead 
to boosts of an order of magnitude or more. In particu- 
lar regions of the parameter space, the dark matter self- 
interaction can drive a model that would be undetectable 
or ruled out by direct searches in terms of er^ 1 to be de- 
tectable at contemporary or future high-energy neutrino 
telescopes. This is an interesting possibility, because this 
implies that the indirect neutrino signal from the Sun 
would not be related to either direct search results or 
indirect neutrino signals form the Earth in a straight- 
forward manner. Each contour of [3 in Fig. [3J exhibits 
a distinct break as it nears the equilibration boundary. 
This is because solutions that are well away from equi- 
librium tend to lie in the linear portion of N x (t), prior to 
any significant opportunity for exponential growth (see 
Fig.[T]). As a result, the insight gained from Eq. (fT4|) fails 
at low cross-sections and large enhancement factors re- 
quire significantly smaller cp 51 at fixed M x than one would 
estimate from the equilibrium assumption. 

Of course, it is likely that the effect of flux enhance- 
ment due to self-capture is negligible; however, it is useful 
to know just how large this effect could possibly be. It 
is simple to make such an estimate and contours of the 
maximum possible flux enhancement /3 max , are shown in 
Figure |4j The "maximum possible flux enhancement," 
depends upon (<7a v ) an d I compute it as follows. At 
each value of M x , I choose the largest value of <r xx that 
is not already excluded by considerations of large-scale 
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FIG. 3: Factors of flux enhancement in the M x -CTp plane. Each panel shows contours of constant relative flux enhancement 
in models of self-interacting dark matter. The flux enhancement contours are the red lines with negative slope labeled by the 
flux enhancement factor j3. Every panel is labeled according to the assumed values of <r X x and (uav) in each calculation. I 
show for reference on the background in each plot contours of constant muon flux at a detector on Earth in the case where 
annihilation happens through W + W~ (solid, with cut-off at the W mass) and the case where annihilation happens through 
t + t~ (dashed). These reference flux levels are computed assuming that there is no significant self-capture of dark matter, 
C s = 0. Shaded regions at the lower ends of these plots correspond to models that are not yet at their equilibrium levels for 
a Sun of age tq = 5 x 10 9 yr. Roughly speaking, current direct dark matter searches constrain the spin-independent dark 
matter-proton cross section to slightly better than up 51 > 10~ 43 cm 2 at about 10 2 GeV [Hoi. l5ll]. 



structure (see the contours in Fig. [5]). I then compute 
the flux enhancements at each point in the Mx-Cp 1 plane 
for a fixed (ctaw) ((caw) = 1CP 27 cm 3 s _1 in this case). 
The enhancement scales approximately as ~ (a^v)^ 1 as 
given in Eq. (fP5|) . 

Fig. [4] already illustrates that some extreme parame- 
ter combinations may be ruled out with contemporary or 
near future limits from neutrino telescopes. Notice also 
that the contours of constant /3 max are very flat functions 



of M x . This is because existing limits on dark matter 
particle self-interactions scale as <r xx cx M x . This is im- 
portant for comparison with direct detection experiments 
which aim to achieve limits on the dark matter-proton 
scattering cross section on the order of cr^ 1 ~ 10~ 44 cm 2 
in the near future [5(| [EH, [EU, [53[ . Absent dark matter 
self-capture, such a limit would indicate that there should 
be no observable high-energy neutrino flux from the Sun, 
but self-capture can clearly modify this conclusion. 
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FIG. 4: Contours of maximum possible flux enhancements 
given existing bounds on dark matter elastic scattering cross 
sections. This figure is similar to the figure panels of Fig. [3] 
However, in this figure, I show contours of /3 max , the maxi- 
mum possible flux enhancement at each point in the M K — (Jp 1 
plane. I compute this maximal boost at each point by setting 
<7 XX to the maximum allowed value at each value of M x (see 
the limits in Fig. [2}. The quantity /3 max is also a function 
of annihilation cross section and this panel shows f3 m&K for 
{<tav} = 10 -27 cm 3 s _1 . From this figure, it is already clear 
that certain combinations of parameters may be ruled out 
with contemporary or forthcoming neutrino telescope data 
(e.g., Ref. Ei). 



V. SUMMARY AND DISCUSSION 

In this paper I have reconsidered the indirect high- 
energy neutrino signal from within the Sun in models in 
which the dark matter particles have significant self in- 
teractions. The influence of self-interactions is that they 
may allow dark matter particles within the Galactic halo 
to be captured within the Sun by scattering off of dark 
matter particles that have already been captured by scat- 
tering off of nuclei within the Sun. For sufficiently large 
dark matter self-interaction cross sections, this can lead 
to a period during which the rate of capture of dark mat- 
ter particles by the Sun grows in proportion to the num- 
ber of dark matter particles already captured by the Sun. 
The number of dark matter particles within the Sun then 
grows exponentially until this increase is stopped by ef- 
ficient annihilation. The net result is that the Sun may 
contain significantly more dark matter than in models 
with no dark matter self interaction and high-energy neu- 
trino signals due to annihilation of these particles may be 
significantly higher as a result. 

In § IIV1 I showed that mild enhancements of a few 
tens to one hundred percent are possible over a wide 



range of viable parameter space that may be probed 
with contemporary and near-future neutrino telescopes 
[H [EE HI [?1 [Zl, III . Significantly larger flux enhance- 
ments of up to a factor of ~ 10 2 are possible in more 
extreme corners of the dark matter parameter space and 
at flux levels that are not yet within reach of near-term 
neutrino telescopes. Ten percent enhancements are not 
particularly interesting at present because intrinsic er- 
rors in the flux predictions are at the tens of percent 
level [IE EH, HE [IE HI] , but this situation may improve 
as experimental advancements drive renewed interest in 
this signal. 

Large flux enhancements require large dark matter self- 
interaction cross sections a xx > 10~ 24 cm 2 and rela- 
tively small dark matter mutual annihilation cross sec- 
tions (<tav) 10 -27 cm 3 s _1 . The small annihilation 
cross section allows the number of dark matter parti- 
cles within the Sun to grow exponentially for a pro- 
longed period of time, which enables large flux enhance- 
ments. More specifically, the neutrino flux enhancement 
grows approximately as cx (7 2 x /<7p I (c7Av), neglecting the 
possibility that for some values of these parameters the 
flux may not reach its equilibrium level for a sun of age 
tq = 5 x 10 9 yr. As a consequence, large enhancements 
require a disparity between scattering and annihilation 
cross sections that may be unfamiliar. However, such a 
disparity is practicable and, in fact, previous proposals 
of self-interacting dark matter rely on just such relative 
differences in cross sections in order to produce signifi- 
cant astrophysical effects without annihilatin g a ll of the 
dark matter in the early universe (e.g., Refs. |22l. l24l|). 

The high-energy neutrino flux enhancement I compute 
may have several interesting implications. In § Mil and 
in the Appendix, I show that the flux from within the 
Earth will not be enhanced due to dark matter self inter- 
actions. In the standard picture, the flux from within the 
Earth can be predicted relative to the solar flux. In vi- 
able contemporary models, the Earth signal is often not 
yet equilibrated and the Sun-to-Earth flux ratio depends 
upon the dark matter particle mass as well as the capture 
and annihilation rates. In the self-interacting scenario, 
the relation between the Sun and Earth neutrino fluxes 
may no longer hold and significant deviations from any 
predicted ratio may be a sign of dark matter self inter- 
actions. 

Likewise, experiments that undertake direct dark mat- 
ter searches may exploit indirect detection methods to 
cross-check limits and/or detections. The correspon- 
dence between direct detection experiments and high- 
energy neutrinos from the Sun is relatively straightfor- 
ward. Though direct detection rates and high-energy 
neutrino fluxes depend on somewhat different integrals 
over the dark matter velocity distribution, in the stan- 
dard picture they both grow in proportion to the product 
of the local dark matter density multiplied by the dark 
matter- nucleon cross section, cx p x ON ■ If dark matter ex- 
hibits considerable self-interaction, this correspondence 
is also broken. The neutrino flux from the Sun may be 
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significantly larger than would be predicted based on the 
limits or detections from direct detection experiments. 
One extreme possibility is that neutrino fluxes that may 
seemingly be ruled out by direct searches (based upon 
limits on cjp 1 and/or (Tp D ) may be realized due to the en- 
hancement from dark matter self-interactions. A broken 
correspondence between the neutrino fluxes and direct 
search results may signal dark matter interactions. 

The pace of the quest to identify the dark matter is 
picking up rapidly. Neutrino telescopes play an impor- 
tant role in this endeavor and the indirect limits from ex- 
isting facilities are already competitive with direct search 
techniques. The indirect, high-energy neutrino signal 
from the Sun may serve as a unique probe of new physics 
confined to the dark sector, and experimental advance- 
ments in the near future should shed new light on the 
properties of the dark matter. 



APPENDIX: CAPTURE AND SELF-CAPTURE 
OF DARK MATTER PARTICLES IN THE SUN 

In the interest of completeness, I give a brief discussion 
of the capture of self-interacting dark matter particles 
within the Sun in this appendix. The treatment here is 
not original, save for the fact that I consider dark matter 
particles interacting among themselves, and follows the 
lucid discussion given in the series of papers by A. Gould 
[Tol . [Til . Il2| . I conclude this section with the rate of dark 
matter particle self-capture. For the results in the main 
text, I use the full formulae of Ref. [l(| to compute dark 
matter particle capture off of nucleons. 

Gould begins by considering capture in an individual 
spherical shell of the body on which capture is occurring 
(the Sun in this case) of radius r and local escape speed 
Vescir). About this shell, consider a bounding surface of 
radius R so large that the gravitational field due to the 
Sun is negligible at R. Let the one-dimensional speed 
distribution function of dark matter particles at this shell 
be /(u), where u is the speed at infinity and the integral 
of f(u) over all speeds gives the number density of dark 
matter particles. The inward flux of particles of speed 
u at angle 9 relative to radial across the surface at R is 
then 



1 



dud cos 2 (0) 4 



f(u)u. 



(A.l) 



Changing variables from cos 2 (6*) to the specific angular 
momentum J = Ru sin(0) , and integrating over the sur- 
face area of the sphere at R gives the rate at which dark 
matter particles enter the surface per unit time, per unit 
speed, per unit angular momentum, 



du dJ 2 



nf(u) 



(A.2) 



Notice that I have written this so that the quantity 71^ 
has dimensions of inverse time. 



Take il(w) to be the rate at which a particle with speed 
w at the shell at r scatters to a speed less than v csc {r). 
Infalling dark matter particles with speed at R of u that 
reach the shell at r, do so with speed 



w = ^u 2 + v 2 sc (r). (A.3) 
The probability of such a particle to be captured is 

dp= OM 2dr Q(m _ J); (A4) 



w 



1 



J 2 



where the quantity 2dr multiplied by the term under the 
radical is the path length through the shell, dividing by w 
converts this to the time spent in the shell, 0(x) is a step 
function and the particular step function above enforces 
the condition that only particles with J < rw intersect 
the shell. Multiplying the rate of incoming particles in 
Eq. (|A.2|) with Eq. (IA.4j) , the differential rate of capture 
within the shell is 



dC 



dr du d J 2 



du dJ 2 
2tt/(u) 



dP 
dr 

n(w) 



(A.5) 
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The integral over J 2 can be performed leaving the cap- 
ture rate per unit speed at infinity, per unit shell volume 



dC f(u) 



dudV 



wCl(w), 



(A.6) 



where I have replaced Anr 2 dr with d^. This gives the 
rate per unit shell volume as an integral over the speed 
distribution at infinity, 



dC 

av 



/(«) 



wfl(w) du, 



(A.7) 



and the task remains to determine Q,(w), perform the 
integration over speeds in Eq. (|A.7[) . and integrate over 
the volume of the Sun. 

The rate of scattering in the shell is simply naw, with 
a the scattering cross section and n the number density 
of targets. The case of most practical interest is velocity- 
independent and nearly isotropic scattering of infalling 
dark matter particles against targets that are effectively 
at rest with respect to the capturing body. In this case, 
the fractional loss of kinetic energy in a given scattering 
event is a uniform distribution over the interval 



AE 4M x m 
< < - 

~ E ~ (A/ x + m) 2 ' 



(A. 



where M x is the mass of the dark matter particle and m 
is the mass of the particle it scatters off of. The dark 
matter particle must lose a fraction of its kinetic energy 
AE/E > u 2 /w 2 in order to be captured. If the condition 



4M K m/(M K + m) 2 > u 2 /w 2 



(A.9) 
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holds, the probability that an individual scattering event 
leads to capture is 



Pcap 



1 - 



u 2 (M x - raf 
v 2 sc (r) 4:M x m 



(A.10) 



Therefore, if Eq. (IA.9|) holds, 



Q(w) = nav csc (r)- 



:(r) 



1 



(M x 



(A.ll) 

At least one property of Eq. (jA.lljl is familiar. Capture 
is most efficient when both projectile and target are of 
the same mass and becomes less efficient as the masses 
become mismatched. 

Combining Eq. (|A.7|) with Eq. (|A.11|) yields the cap- 
ture rate per shell volume in the Sun, 



dC 
dV 



nav„ 



/(«) 



(M x 



du. 



(A.12) 

Gould has evaluated this expression for the case of a 
Maxwcll-Boltzmann speed distribution including possi- 
ble form-factor suppression of scattering with large nuclei 
at high momentum transfer [ToL [Tl|. However, the gen- 
eral formulae are rather unwieldy, the integrations are 
lengthy but straightforward, and presenting them does 
not add significantly to the insight needed for my pur- 
poses. As a result, I will not present the general formulae 
and will move to a particularly simple special case. 

Of particular interest for the present paper is the cap- 
ture of dark matter particles in the halo by other dark 
matter particles that have already been captured within 
the Sun. As a consequence, I will evaluate Eq. (|A.12[) for 
the special case of m = M x and for capture by the Sun 
moving with speed Vq = 220 kms -1 through a Maxwell- 
Boltzmann distribution of dark matter particles with dis- 
persion v — 270 kms - . The distribution function can 
then be written 



fix) 



2n x 2 _ x 2 a sinh(2a;?7) 
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(A.13) 



in terms of the dimensionless variables x 2 = 3(u/v) 2 /2 
and T) 2 = 3(vq/v) 2 /2. Integrating over the speed distri- 
bution yields 



dC 
dV 
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2 V T] 



(A.14) 



The total capture rate now requires integrating over 
the volume of the Sun. This gives 
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^csc (R Q ) erf (77) 
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(A.15) 



The last integral can be re- written conveniently by defin- 



ing a dimensionless potential 



;(r)A4c( R 0). in 



which case the last integral is the product of the total 
number of targets N and the average of <p over all tar- 
gets within the Sun, 



n /3 (vt v WoscCHq) ,7, erf(?7) 

C = \ n nx<7v csc (RQ) N{(p) . 

V 2 v i~i 



(A.16) 



The numerical factor of 1/3/2 in Eq. (|A.16[) differs from 

the factor \/6/ir given in Refs. [l^, HH, G3 because Gould 
defined the error function erf(cc) with an unconventional 
normalization. 

I have assumed that M x = m to derive Eq. (|A.16[) . 
However, so long as the mass of the target and and pro- 
jectile are not very mismatched, scattering will be likely 
to lead to capture and Eq. (|A.16|) will be a relatively good 
approximation for the capture rate. In the case of cap- 
ture by scattering off of nuclei, the relevant cross section 
is the elastic scattering cross section off of the nucleus of 
interest a = ctn and N is the number of such nuclei in 
the Sun. The total capture rate due to scattering off of 
all nuclei is the sum of the individual rates for all of the 
different nuclear species within the Sun. 

For dark matter self-capture, the relevant cross sec- 
tion is the elastic scattering cross section of dark matter 
particles with themselves a = <7 XX and N — N x is the 
number of dark matter particles already captured within 
the Sun. Therefore, the dark matter self-capture rate co- 
efficient referred to in the main text can be approximated 
as 



,-. /3 , ^ CS c(R Q ) ,7 v erf (77) 
U, = \l- n x cr xx v osc (R j = {<p x ) . 

2 V T] 



(A.17) 



As discussed in the text, captured dark matter particles 
typically occupy a very small range of radii within the 
Sun (typically confined to only a few percent of R©), in 
which case (0 X ) ~ 5.1 [Toj . 

In the case of dark matter particle self-capture, there 
is one additional complication that must be accounted 
for that is not relevant for capture off of nuclei. The 
Sun is optically-thin to the propagation of dark matter 
particles, so a target dark matter particle that receives 
too much kinetic energy relative to the solar core will be 
ejected resulting in no net gain of dark matter particles. 
Therefore, not only must the collision result in an energy 
exchange of AE/E > u 2 /w 2 , but it must be limited to 
AE/E < v 2 sc (r)/w 2 . This modifies the capture proba- 
bility per collision (again, taking m = M x ) to 



Pcap — I 

V 

and the capture rate to 

n(to) 



e(v csc (r)-u) (A.18) 
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-(v 2 esc (r)-u 2 ). 



(A.19) 
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When v csc (r) 3> u, this modification is relatively minor. 
This is because in this situation, the incoming dark mat- 
ter particle must only lose a small fraction of its total 
energy to be captured and does not necessarily impart 
enough energy to escape on the target dark matter parti- 
cle. This is generally the case for the Sun, because escape 
from the solar interior requires speeds at least two times 
larger than the typical speed at infinity of a dark matter 
particle. However, the escape speed from the Earth is 
significantly smaller than the typical speeds of dark mat- 



ter particles, so collisions within the Earth that lead to 
capture of the infalling particle will almost always lead 
to ejection of the target. In fact, most interactions of 
this kind will lead to both infalling particle and target 
being unbound from the Earth. An interesting question 
is to ask whether self-interactions may scour the Earth of 
captured dark matter particles, but a comparison of the 
relevant rates along the lines leading to Eq. (JT5J) in § IIIII 
shows that the removal rate is significantly less than the 
capture rate for parameters of interest. 



J 



This small modification results in a significantly more complex formula for the rate of capture. The calculation 
follows according to the simple estimate given above. Again, the integrations are lengthy but straightforward, so I 
will only quote the result. The full rate of capture accounting for the potential recoil and ejection of the target dark 
matter particles is 



- n x cr xx Ve SC (R©j = V 

2 v 




((^ x erf (x v + rj)) - (4> x erf{x v - 77))) 



[(erf(x v + r))) - (erf(x v - 77)}] 



+2e"" 2 - [(e(^ 2 ) - (e^+^ 2 )] + -J(0 )V ) - -(J(x v - ^x Y + V )) 



(A.20) 



where x\ = 3(v esc (r)/v) 2 /2, the brackets about a quantity, such as "(9)," designate the average over all captured 
dark matter particles of the quantity q, and the integral J(s : t) = f q 2 e~ q dq. The first term in this relation is the 
simple result from Eq. (|A.17[) . The second term in the first set of square braces results from truncating the integral 
over the speed distribution at u = v csc (r). Typically, x v ± r) > 1, so this term will be small in comparison to the first 
term. The terms within the curly braces come from the new piece in the capture rate Eq. (|A.19[) . The factor that 
multiplies the terms in curly braces is typically of order ~ 0.06 for the Sun. Consequently, the new terms in Eq (|A.20|) 
collectively represent relatively small modifications to Eq. (|A.17|) . This fits the heuristic understanding that ejection 
due to recoil will be important only when v QSC (r) < v ~ v@. 

I 



Though the above sketch of Gould's derivations is in- 
structive for present purposes, the formulae I present here 
do not suffice to make an adequate estimate of capture 
by nuclei within the Sun. In all of the detailed results 
in § IIV1 I use the full formulae given in Ref. [l0| and re- 
peated in the review of Ref. [l[. I tak e = 220 kms -1 , 
v = 270 kms -1 , p x = 0.4 GeV/cm 3 '77j], the solar mass 
distribution of Ref. [l(|, and the elemental abundances 
given in the review of Ref. [78[ . I show a specific example 
of my calculations of the rate of capture of dark matter 
particles from spin-independent scattering off of nuclei 
in the Sun with a spin-independent cross section for dark 
matter-proton scattering of <7p = 10~ 43 cm 2 in Fig- 
ured] In addition to the total capture rate, I show also 
in Fig.[5]contributions to the total capture rate from scat- 
tering off of several of the most important nuclei within 
the Sun. Capture off of Hydrogen is down by roughly 
two orders of magnitude throughout most of this range 



due to the lower cross section relative to heavier nuclei 
and the unfavorable scattering kinematics for heavy dark 
matter particles. 
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M, [GeV] 



FIG. 5: Capture rates of weakly-interacting dark matter par- 
ticles used in the calculations in the main text. In this panel, I 
assume ap 51 = 10 -43 cm 2 and I show capture rates as a func- 
tion of dark matter particle mass. In addition to the total 
capture rate, Cj OTAL , I also show capture rates off of several 
elements within the Sun for those elements most important to 
capture via a scalar interaction. These are He (C He ), the sum 
of C, O, and N (C CNO , Oxygen is the most important of the 
CNO elements individually), the sum of Fe and Ni (C Fc+Ni ), 
and Ne (C Nc ). At high-mass, the capture rate approaches 
C c oc M^ 2 as expected [ljj. 
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